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Spectral Energy Distributions of Gamma Ray Bursts 
Energized by External Shocks 

Charles D. Dermer^, Markus Bottcher^'"^, & James Chiang'* 
ABSTRACT 



Sari, Piran, and Narayan have derived analytic formulas to model the spectra from 



o 

^-v ■ gamma-ray burst blast waves that are energized by sweeping up material from the surrounding 

^S| , medium. We extend these expressions to apply to general radiative regimes and to include 

^H ' the effects of synchrotron self-absorption. Electron energy losses due to the synchrotron 

■ self-Compton process are also treated in a very approximate way. The calculated spectra are 

compared with detailed numerical simulation results. We find that the spectral and temporal 
1/^ ' breaks from the detailed numerical simulation are much smoother than the analytic formulas 

imply, and that the discrepancies between the analytic and numerical results are greatest 
S^ , near the breaks and endpoints of the synchrotron spectra. The expressions are most accurate 

^SJ ' (within a factor of ~ 3) in the optical/X-ray regime during the afterglow phase, and are more 

r^ ■ accurate when Ce, the fraction of swept-up particle energy that is transferred to the electrons, 

is ;S 0.1. The analytic results provide at best order-of-magnitude accuracy in the self-absorbed 

radio/infrared regime, and give poor fits to the self-Compton spectra due to complications from 
0^ ■ Klein- Nishina effects and photon-photon opacity. 

ON 



Introduction 



Oh' 

o 

-*— > ■ The discovery of X-ray and long wavelength afterglow radiation from gamma ray bursts (GRBs) finds 

j^ I a convincing explanation through the blast-wave model (e.g. Paczyhski & Rhoads 1993; Meszaros & Rees 

1997). In this model, a coalescence or collapse event, possibly involving the formation of a black hole, 
produces a weakly baryon-loaded fireball with directional energy releases dE/dfl as large as ~ 10^^ ergs 
r> ' sr~^. The blast wave sweeps up material from the surrounding medium during its expansion. The swept-up 

C^ ' particles are a source of free energy in the comoving blast wave frame. Depending on the magnetic field 

strength and the intensity of the surrounding radiation field, this energy is primarily radiated through 
synchrotron and Compton processes. Particle escape, or other channels of particle energy loss such as 
adiabatic losses or secondary production for the high-energy protons (Pohl & Schlickeiser 2000), may 
dominate at certain times in the blast wave evolution and in certain regimes of particle energy. 

Sari, Piran, & Narayan (1998; hereafter SPN98) have derived useful formulas to model the emissions 
from a decelerating blast wave in the afterglow phase, and asymptotic forms were given in the adiabatic and 
radiative limits of blast-wave evolution. These formulas have been widely used to interpret observations of 
GRB afterglows. Their work was restricted to the optical and X-ray range where synchrotron self-absorption 
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(SSA) is generally unimportant. Moreover, the synchrotron sclf-Compton (SSC) process was not treated, 
and no comparison of the analytic expressions with detailed numerical simulations was carried out to assess, 
for example, the inaccuracy incurred by not considering adiabatic losses on the evolution of the electron 
energy spectra. 

In this article we generalize the formulas of SPN98 to include SSA and SSC processes and to apply 
to arbitrary radiative regimes in the evolution of the the GRB blast wave. The results are compared with 
a detailed numerical model (Chiang & Dermer 1999; Dermer, Chiang, & Mitman 1999). The numerical 
treatment approximates the blast wave as a thin shell within which the energized particles radiate in a 
uniform magnetic and photon field. It therefore applies to the forward shock, which is the appropriate 
system that should be compared with the analytical formulas. Emissions from the reverse shock are not 
considered here. 

In the next section, we summarize the parameters that enter into the standard relativistic blast-wave 
model for GRBs, and describe the electron distribution function and evolution of the blast wave. In Section 
3, the SSA process is considered. Section 4 describes our treatment of SSC processes. A presentation of the 
results of the semi-analytic model, including a comparison with numerical data, is provided in Section 5. 
We summarize in Section 6. 



2. Electron Distribution and Blast Wave Evolution 

We adopt the notation of SPN98, except where indicated. The power injected in nonthermal electron 
energy in the comoving blast-wave frame is assumed to be equal to a constant fraction e^ of the power 
E of particles swept-up by the blast-wave. For an uncoUimated blast wave expanding isotropically into 
a uniform medium with density n, E = 47rr^nmpC'^/3(r^ — F) (Blandford & McKee 1976). Here F is 
the blast-wave Lorentz factor and cf3 — c(l — 1/F^)^/^ is its speed. The electrons are assumed to be 
injected in the form of a power law with index p and with minimum and maximum Lorentz factors 7™ and 
72 = Cmax • 4 X 10^/[i?(G)]^/^, where the latter expression derives from a comparison between the maximum 
electron acceleration rate and the synchrotron loss rate, and emax < 1- We set Cmax = 1 in the examples 
shown in this paper. If the number of energized nonthermal electrons is equal to the number of swept-up 
electrons, then 

^ CtZ~\1 ) = 1 + ^e(r - l)m,/m, , (1) 

P - ^ 7m - 72 

which reduces to the commonly used expression 7™ = edip ~ 2)/(p — l)](r — l)mp/me when 72 3> 7m 3> 1 
and p > 2. Eq. (nl) is easily solved numerically to give an accurate value for 'jm- Note that the meaning of 
Ee in the numerical simulation differs from its meaning in the analytic model (Wijers & Galama 1999); in 
the latter case, e^ represents the ratio of (residual) nonthermal energy in electrons to that in protons. 

Following the reasoning of SPN98, electrons will cool in the comoving fluid frame by synchrotron and 
Compton losses to Lorentz factor 

4XmpnceB(Jt'r^t{l + /KN^ph/ws) 

by observing time t. Here A is the compression ratio (set equal to 4 in the following calculations), and 
ub = B'^/Stt = XnmpC^eB^^ and Uph are the magnetic field and photon energy densities, respectively. 
When Klein-Nishina effects on Compton scattering are important, Mph is corrected for the reduction in the 
scattering efficiency using the factor /kn (< 1), as discussed in more detail below. 



The distribution of energized electrons is approximated by 



s-i J 7 ^ , for 7o < 7 < 7i 



where Nf, — Anr^n/^ is the total number of nonthermal electrons. In the slow cooling limit, 70 = jm, 
7i = 7c, and s — p, whereas in the fast cooling limit 70 = 7c, 71 = 7mi ^^nd s = 2. The blast wave's radiative 
efficiency e = eeEradj where e^ad = 1 in the fast cooling limit, and Ej-ad = (7m/7c)^~^ in the slow cooling 
limit (Moderski, Sikora, & Buhk 1999). 

The evolution of T in radiative regimes bridging the adiabatic and radiative limits can be quickly 
obtained by numerically solving the equation of motion dV /dm — — (F^ — 1)/M . Here dm = Anr'^mpudr 
and M — J dm [e + r(l — e)] is the total energy in the blast wave, including both rest mass and kinetic 
energy. This prescription for the change in M is consistent with energy conservation at all energies, and 
reduces in the nonrelativistic limit to Oort's solution based on momentum conservation rather than the 
Sedov self-similar solution, as noted by Huang, Dai, & Lu (1999). Analytic solutions for these equations 
when e is constant are presented by Bottcher & Dermer (2000), where observational implications of general 
radiative regimes are discussed. 



3. Synchrotron Self- Absorption 

In the absence of SSA, the analytic approximation to the synchrotron radiation consists of a low-energy 
hard power-law emission spectrum with spectral luminosity L^, oc v^'^ at v < Vq, connected to an 
intermediate branch with L^ oc j/^('*^i)/2 at j/q < i^ < t^i, which joins to a higher energy branch with 
L,y ex i/~P" at 1^1 < 1/ < 1^2- The frequencies of the breaks in the synchrotron spectrum are given by 
Vi = rjfeB/{2Trmec), with i = 0, 1,2. (Observed frequencies should be divided by the factor (1 -I- z) to 
implement a redshift correction for a sourceat redshift z.) When SSA is important, the spectrum exhibits 
a break at 1/ < Va, where Va is the self-absorption frequency. If I'a < fo, then F^ (x ly^ a,t v <Si Va- When 
fa > ly'o, Fij (X ly^ at frequencies ly < vq and F,y ex i^^" for fQ < ly < i/a- We do not treat the case h'a > t^i, 
though it is a straightforward generalization of the results presented here. 

When J/ < fo, we calculate the SSA coefficient from the standard expression 



djP{,.',j)^' — [^^], (4) 

where ^^(7) = (s — l)Afe7o "'^7 ** from eq. (H), the blast-wave volume Vbw — 47rr^A', and A' is the 
comoving- frame width of the shell. In the specified limit, the calculation of the absorption coefhcient is 
simplified because we can use the asymptotic form for the synchrotron emissivity function 

S^/^e^Bsina^, , , Aire^B . v .w, , , 

where the right-hand-side applies in the limit v <^ v^r, where v^r — Zvb^^ /'i- is the critical frequency. The 
exact expression for F{x) is given by, e.g., Rybicki & Lightman (1979). Note that r(l/3) = 2.679 is a 
Gamma function. In this expression vb = e_Bsina/27r?7ieC and we let the pitch angle a — 7r/2. Substituting 
eq.(H) into eq.(0)and letting the relation k^i/^' = 1 define the comoving- frame absorption frequency, we 
obtain for the observed SSA frequency the expression 

fa(Hz) = 9.9 [ii±^fc-il]3/5 (,n)3/5s2/5 1^ ^^^^ ^^^^^ (g) 

(s + 116) 7o 



All of the terms in eq.(^ are in cgs units. 

For the case I'o < lya < t^i, we use the results quoted by Gould (1979) for the absorption coefficient. 
Thus 

K,, = c{s)rlK{uJu'){vBli^'t^^'^'^ , (7) 

where Uq — c/rg and re is the classical electron radius. The function c{s) is tabulated in Gould (1979), and 
fee = (s — l)-/Ve7o~^/''4w — {s — 1)7-7170"^ /(3A'), using eq.(^). Defining the absorption frequency in this 
case by the value of v satisfying dL„/dv = 0, noting that L^y oc v^^'^\V — exp(— i)] and t = Kjy' A', we obtain 
for the observed SSA frequency the expression 



^ (■s-l)c(s)rgm7g Vo ^2/(s+4) ,/j<+2)/(s+4) 
3t. 



z/^(Hz) = \- ^^ ]^/is+''-i jy^^ ii\ ^ I Y ,when vq <Va <v\ . (8) 



The term i^ is the solution to the transcendental equation exp(— t^) = 5/[5 + (s + 4)is] derived in the slab 
approximation. Representative values are t^ = 0.35,0.50,0.61 for s — l^hjl^ and 3, respectively. 

The synchrotron spectral luminosity -Lzy.max at frequency vq is given by ii^^max = N efne,(? ot^ B j (ie) 
(SPN98). This holds irrespective of the level of Compton losses, because I'o^i^.max represents the measured 
synchrotron power from A^e electrons with Lorentz factor 7 ~ 7o- Consequently, the synchrotron spectral 
power when Vg, < i'q is 



for Va <V < Vq, 

for i^Q < u < vi 




E_\-{s-i)/2 f„„ ,. ^ ., ^ .. ' y^J 



and Pa is given by eq. (g|) . When vq < Va < 1^1, Va is given by eq. (g) , and the first two cases in eq. (^ should 
be replaced by 

r(£0)5/2(i^)-(s-l)/2(Ji)2 fori.<j.o 

The last two cases of eq. ^ also hold at frequencies v > Va when vq < Va < vi- 



4. Synchrotron Self-Compton Emission 

Moderski et al. (1999) show that Compton processes are energetically negligible when e = egErad ^ Ae^ 
(see also Sari, Narayan, & Piran 1996). This is a conservative estimate because it assumes that all scattering 
takes place in the Thomson regime. Consequently SSC processes can only be neglected with certainty if the 
energy transfer from protons to electrons is small and the magnetic field energy density is near equipartition. 
To explain the prompt gamma-ray burst emission in an external shock scenario, it is necessary to have 
es ^ 1 in order to be in the weakly cooled regime, in accord with measurements of the hard X-ray/soft 
gamma-ray spectral indices in GRB spectra during the prompt phase (see discussion in Chiang & Dermer 
1999). SSC emission must therefore be treated in this case, and also when eg ^ 0.1, which is an important 
regime to consider for good radiative efficiency in GRB blast waves. 

We employ a very simple treatment for the SSC process by approximating the synchrotron emission 
by a (5-function distribution at the comoving frame frequency vi/T. This is the frequency of the peak of 
the power output if 2 < p < 3. As noted earlier, the relative energy radiated in Compton and synchrotron 



photons is, neglecting Klein-Nishina effects, equal to Uph/uB — Lsyn/('Us47rr^c), where Lgyn is the 
synchrotron power in the comoving frame. The synchrotron power is 



72 

2 



4 r^ 

Lsyn = o CCTtUb / d^N{'^)^ 



3 3 — s2— p 2— p 3 — s 

Because of the strong reduction in scattering efficiency when photons are scattered in the Klein-Nishina 
limit, we shut off the integration when ^hvi/ {rme(?) > 1, or when j — 'j — TnieC^ /{hiyi). This implies that 
no photons with frequencies i^kn > (TmeC^ /h)'^i'^ are emitted. The energy density of photons that are 
efficiently scattered in the Thomson regime is obtained by repeating the exercise in eq. (O) , but integrating 
over the range 7o < 7 < 7 rather than over 70 < 7 < 72- This yields a ratio /kn of the Compton power 
which approximately takes into account Klein-Nishina effects to that without Klein-Nishina effects. This 
ratio is roughly given by /kn = 1 if 7/71 > 1, and /kn ~ (7/71)"^ '' if 7/71 < 1- Because 7 is the electron 
Lorentz factor for which Klein-Nishina effects become important, it is also necessary that 7 > 70 for 
significant Compton emission. 

The spectral power of the Compton-scattered radiation at frequency i'^ — 791^1 is 

3 V^ Ub Vq Ub 

The SSC component in this approximation is therefore given by 

r {^f'\ for V < v^ 

^? = iS„.a. X (|)-(^-^)/^ for .0^ < . < .f ^ 7^1 , (13) 

[ (^)(^"'Hi^)-^/' for v?<v< v^ = ^Iv, 

provided i^kn > ^2- When this is not the case, the SSC emission is set equal to zero at v > i^kn- 



5. Results 

Figure 1 compares spectral energy distributions (SEDs) obtained from the semi-analytic formulas 
derived above with those calculated using a detailed numerical simulation of GRB blast waves. The 
numerical model is described by Chiang & Dermer (1999) and Dermer et al. (1999). The parameters in 
this calculation represent a standard set chosen to produce a bright, prompt GRB with a peak gamma-ray 
luminosity ^ 10^^ ergs s~^, a duration ^ 10 s, a vF^, peak synchrotron frequency i^pk corresponding to 
photon energies of several hundred keV, and a hard spectrum approaching a = —1/3 at i^ ^ fc'pk, where 
the flux density F{v) oc i/~". For this calculation, a fireball with total energy Eq — 10^'* ergs expands to 
produce an uncollimated blast wave (effects of blast-wave collimation are numerically treated by Moderski 
et al. 1999, Panaitescu & Meszaros 1999, and Dermer et al. 1999). For the reasons given above, a very 
small magnetic-field parameter of e^ — 10^* is used. The injection index p — 2.5 is used to match spectral 
indices inferred from afterglow spectra at optical and X-ray frequencies. The curves are labeled by the base 
10 logarithm of the observer time, and the numerical results are highlighted by data points. 

As can be seen from Fig. 1, the behavior of the blast wave emissions at soft gamma-ray energies and 
below is qualitatively reproduced by the analytic expressions. The synchrotron emission reaches a peak at 



early times and at high photon energies. The deceleration of the blast wave causes the observed synchrotron 
fluxes and i^pk values to decrease with time. The fit is best during the afterglow phase at f 0^ s ^ f ;S f 0^ s 
at frequencies 10^'* Hz ^ v ^ 10^^ Hz. The characterization of the high-energy gamma-ray SSC component 
is very poor, especially during the prompt phase. This is due to the oversimplified way of dealing with 
the Klein-Nishina reduction in the cross section. Progress in improving the analytic representation of the 
SSC emission comes at the expense of cumbersome formulas and which furthermore still omit the effects of 
photon-photon attenuation that can be important at TcV energies. Thus we do not go beyond the simple 
analytic treatment of the SSC component given in the previous section. 

The analytic SEDs also provide a poor representation of the hard X-ray spectra shortly after the 
prompt emission phase at i ~ 10-10'^ s. This discrepancy has at least two causes. The first is the neglect of 
adiabatic losses on the electrons. These losses degrade the low-energy portion of the electron distribution 
function and produce a smoother low-energy cutoff in this distribution function than assumed in the 
analytic representation. The second is the analytic approximation for the elementary synchrotron emissivity 
spectrum, which consists of two connected power laws. This shortcoming is illustrated in Fig. 2, where we 
plot synchrotron emissivity spectra using the formulas of Crusius & Schlickeiser (1986), which have been 
averaged over pitch-angle and magnetic field direction, for power-law electron distribution functions with 
index s. These results arc compared with an analytic representation of the synchrotron emisivity spectrum 
derived in a (5-function approximation (Dermer, Sturner, & Schlickeiser 1997). The J- function approximation 
yields the expression of SPN98 for the synchrotron emission except for an overall normalization of a factor 
2/(p — 1). Fig. 2 shows that the analytic results overestimate the flux by as much as a factor of ^ 4 at 
the low-energy break where the analytic power-law connects to a power law with energy index a ~ ^1/3 
extending to lower frequencies. At the high-energy endpoint of the synchrotron spectra. Fig. 2 shows that 
even larger discrepancies between the numerical and analytic results are obtained, though these are usually 
concealed by the SSC component. 

The effects of adiabatic losses on the low-energy distribution function becomes less important at later 
times. Thus the flux discrepancies in the low-energy portion of the synchrotron spectra decrease with time, 
but they still can differ by an ordcr-of-magnitude. 

Fig. 3 shows the effect of considering a stronger blast-wave magnetic field in the model. All parameters 
are the same as in Fig. 1, except that €b = 0.1. With this stronger field, the GRB would produce a brighter 
observed synchrotron emission spectrum, though most of its power would be carried by photons with higher 
energies. In this case, the bulk of the synchrotron power is emitted in the form of photons with energies of 
~ 50 MeV. Note, however, that this model GRB would be dimmer in the BATSE triggering range than the 
model GRB shown in Fig. 1. The strength of the SSC component is much weaker than in the model of Fig. 
1, particularly in the afterglow phase. Although the overall behavior of the numerical simulation results 
are reproduced by the analytic model, large discrepancies still persist for the same reasons as previously 
discussed. In fact, the analytic model at optical and X-ray energies with SSC effects included is not much 
better than a model which neglects them, as indicated by the 10^, 10® and 10^ s curves. 

The analytic representation of the SEDs improves when smaller values of e^ are used. Fig. 4 shows 
the effects of reducing eg to a value of 0.1, with all other parameters the same as before. This has the 
effect of further reducing the relative importance of the SSC component. Nevertheless, order-of-magnitude 
discrepancies between the analytic and numerical models remain at the endpoints of the spectra. 

Figs. 5a and 5b compare analytic and numerical SEDs for a model GRB with parameters derived by 
Wijers & Galama (1999) to fit afterglow data from GRB 970508. Here the total energy injected into the 



fireball is E = 3.5 x 10^^ ergs, which is considerably smaller than the value used in the previous cases, and 
the density of the external medium is very dilute with n = 0.03 cm""^. Consequently the SSC component 
becomes even less important than before. The analytic expressions provide a reasonable representation of 
the numerical model in the optical/X-ray regime, though the discrepancies at radio frequencies remain 
significant, as is clear from Fig. 5b. This model also does not explain the flux of GRB 970508 in its prompt 
phase, which is much brighter than calculated here. 

The light curves from the semi-analytic and the numerical models for the Fig. 5 parameters are shown 
in Fig. 6 in an L^, representation. For this example, there is good agreement between the optical light 
curves. The temporal breaks due to the peak of the synchrotron emission spectrm from the lowest energy 
electrons passing into the specified waveband are much smoother than the analytic model, and discrepant 
in flux by an order-of-magnitude. The cooling break, which is clearly defined in the V band light curve at 
t Ri 10^ s in the analytic model, is almost indistinct in the numerical calculation. 



6. Discussion 

We have examined the accuracy of analytic formulas used to model GRB afterglow radiation by 
comparing results from these expressions with calculations obtained with a numerical simulation model that 
treats the most important processes in relativistic blast waves. The advantage of the analytic model is that 
it can be easily used to model data and extract parameter values, but its accuracy is uncertain. By contrast, 
the numerical model contains all relevant processes and should therefore be used to fit multiwavelength 
GRB data, but the available numerical model (Ghiang & Dermer 1999) requires ~ an hour (for eb ^ 1) to 
a day (for cb ~ 0.1-1) to run on a SparcStation. This is too long to explore parameter space as is required 
when fitting data. For fits to data, it seems most useful to establish general parameters using the analytic 
model, and then adjust parameters in the numerical simulation to obtain a precise fit to data. 

We find that the analytic results are reasonably accurate (within a factor of ^ 3) for modeling 
X-ray/optical data in the afterglow phase (i.e., several hours to many days after the GRB event), but is 
inaccurate in the prompt and early afterglow period, especially when Cg ^ 0.1 and Ce/es 3> 1. Endpoint 
effects in the IR/radio and gamma-ray regimes, and the simplified treatment of the SSG process which 
is important at gamma-ray energies, limit the accuracy of the analytic model to no better than an 
order-of-magnitude. As expected, the spectral and temporal breaks in the analytic approxmations are much 
sharper than found in the numerical calculations. The use of the analytic model to provide simple fits to 
optical and X-ray afterglow data is generally warranted, but the analytic model should be applied with 
caution to radio data or optical and X-ray data in the prompt and early afterglow phases. 

We thank Professor Peter Meszaros for comments on the manuscript. This work is supported through 
the NASA Astrophysical Theory Program (DPR S-13756G) and the Office of Naval Research. The work of 
M.B. is supported by NASA through Ghandra Postdoctoral Fellowship grant number PF9-10007 awarded 
by the Ghandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for NASA 
under contract NAS8-39073. 
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Fig. 1. — Calculations of vLi^ spectral energy distributions emitted by a relativistic blast wave with initial 
Lorentz factor To = 300. Curves are labeled by the base 10 logarithm of the observing time in seconds. 
The analytic and numerical models are shown, with the latter curves identified by data points. The other 
parameters of the calculation are shown in the legend. 
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Fig. 2. — Comparison between a i5- function approximation to the synchrotron emissivity spectrum for a 
spherical emitting region of radius R = 10^^ cm, and the synchrotron emission obtained by averaging over 
pitch angles and magnetic field directions. Curves are labeled by s, the index of the steady-state number 
distribution of electrons, which are assumed to be distributed as a power law between 7niin = 100 and 
7max = 10^. The total energy of the electrons is E — 10^° ergs, and the mean magnetic field strength is 
_B = 1 G. The low-energy power-law branch of the analytic approximation is not shown. 
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Fig. 3. — Same as Fig. 1, but with the magnetic equipartition parameter es = 0.1. Two analytic models, 
one with SSC (solid curves) and one without SSC (dotted curves) processes included, are shown. The curves 
are identical for the t = 1 s case. 



12 



52 



50 



48 



w 

O) 
<D 

-, 46 

o 

o 44 



42 



40 




10 



15 



20 

log^Jv(Hz)] 



25 



30 



Fig. 4. — Same as Fig. 1, but with eb = 0.1 and €e = 0.1. 
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Fig. 5. — Same as Fig. 1, but with parameters taken from the fit to GRB 970508 by Wijers & Galama 
(1999). (a) Comparison of SEDs in a vLi, representation, (b) Comparison of SEDs in an L^ representation. 



14 



32 
31 

^ 30 

'n 
X 

T 29 
w 
w 
P 28 

^ 26 
25 
24 





1 1 1 


1 1 1 1 


-l' 


\>^ \ 


. ...--'S GHz;.._ - 


> 


/prT--'' J^ 




7 


"n % -'*' 


^\. -•' ' 


75 


f^\ ,'-;' 


•■ vK. "^ 


V 




^v 


/ 


/•"\---'\ 


\\Vband - 


/_. 


.-* \ 


^^Ni 




\ 


\N. 




\ 


V\. 


7 


\ ■ 


\\ 




\ 


*.',4 keV \\ 




X 


'• 1 \ \ 


^ 


\ 


\\ " 




130keV « 


'^ \. \y, 










1 1 1 


1 \ ■^- \ '-M 



1 



log^ Jtime (s)] 



Fig. 6. — Analytic and numerical light curves calculated at radio, optical, X-ray and soft gamma-ray energies 
for the case shown in Fig. 5. 



